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INTERINDIVIDUAL VARIATION IN MITOCHONDRIAL ENZYME ACTIVITY IN MALE 
MONARCH BUTTERFLIES, DANAUS PLEXIPPUS L. (NYMPHALIDAE) 
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ABSTRACT. We determined activity levels for the mitochondrial enzyme succinate dehydrogenase in male monarch butterflies overwin¬ 
tering on the central coast of California from October 2001 through March 2002. Mitochondrial activity is important for generating the energy 
required for metabolically demanding activities such as flight, which is essential for monarch reproduction and survival. There is a genetic com¬ 
ponent to variation in flight performance in various species, and individual variation in mitochondrial activity may contribute to these differ¬ 
ences, since mitochondrial enzyme levels often correlate with performance abilities. To understand possible functional consequences of mito¬ 
chondrial activity, it is first necessary to determine the degree of individual variation within the population. We found a high degree of 
interindividual variation in enyzme activity in male monarchs, at least a twelve-fold difference between the lowest and highest activities mea¬ 
sured, with a coefficient of variation of forty-seven percent. In addition, we investigated possible correlations with season, body weight, body 
size, and wing damage. Although there were some month-to-month differences, individual variation in mitochondrial enzyme activity is not ex¬ 
plained by seasonality or body size, and is not related to the degree of wing damage. The results suggest that interindividual differences in mi¬ 
tochondrial enzyme activities are considerable, and worth investigating as a factor in individual performance and success. 
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Monarch butterflies are known for their dramatic 
migration across the U.S.A. and overwintering in mas¬ 
sive aggregations along the California coast and in cen¬ 
tral Mexico (Tuskes & Brower 1978, Brower 1985). 
Flight is energetically expensive; indeed, insect flight 
muscle is notable for the highest metabolic rates and 
power output of any animal tissue (Sacktor 1976, 
Suarez 2000). In addition to long flights of migration 
during fall and spring which may cover thousands of 
miles (Brower 1985), mating involves a short energet¬ 
ically intense nuptial flight. Nuptial flight or mate 
transport to the nearby tree canopy occurs after a male 
successfully couples with a female (Shields & Emmel 
1973). Regulation of metabolism is also important for 
survival, since monarchs rely primarily on stored lipids 
during the overwintering period (Chaplin & Wells 
1982, Masters et al. 1988, Alonso-Mejia et al. 1997). 
Understanding metabolic variation in monarchs could 
therefore provide valuable insight into differences in 
flight performance, reproductive success and survival. 

Relatively little is known about monarch flight mus¬ 
cle metabolism or about intraspecific metabolic varia¬ 
tion in general. Individual variation in flight perfor¬ 
mance has been shown to be partly genetic in moths 
(Parker & Gatehouse 1985). Correlations between 
flight activity or capacity and metabolic differences 
have been seen in Agrotis ipsilon moths (Sappington 
et al. 1995), Colias spp. butterflies (e.g.. Watt et al. 
1983), Epiphyas postvittana moths (Gu 1991) and 
Drosophila melanogaster flies (Barnes & Laurie- 
Ahlberg 1986, Marden 2000). In Colias butterflies, ge¬ 
netic differences in metabolic enzymes also correlated 
with differing mating success by males (Watt et al. 
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1985). There is evidence that monarch butterflies may 
have individual differences in flight ability at cool tem¬ 
peratures (Hughes et al. 1992). Further investigation 
of metabolic parameters may reveal some of the mech¬ 
anisms underlying this individual variation. 

Mitochondrial activity is a critical aspect of metabo¬ 
lism; insect flight muscle metabolism is strongly aero¬ 
bic, so most of the energy is produced by the mito¬ 
chondria (Sacktor 1976). Mitochondrial enzyme levels 
appear to be a good indicator of mitochondrial density 
and overall aerobic metabolism measured by oxygen 
consumption or oxidative capacity in a variety or or¬ 
ganisms (e.g., Holloszy & Booth 1976, Spina et al. 
1996, Putnam & Bennett 1983). The enzymes most 
commonly used as markers of mitochondrial activity 
are succinate dehydrogenase (SDH) and citrate syn¬ 
thase (CS), mitochondrial-specific enzymes with criti¬ 
cal roles in the tricarboxylic acid cycle and therefore 
directly involved in the generation of ATP. 

Individual mitochondrial enzyme activities can vary 
considerably and the variation often correlates with 
performance. The largest interindividual variation re¬ 
ported is an approximately twenty-fold difference seen 
in CS activity in leg muscles of 20 toads ( Bufo mari- 
nus) (Longphre & Gatten 1994). More commonly, dif¬ 
ferences of three- to five-fold are seen in CS activity 
among individuals in a variety of non-insect species 
(e.g., Walsberg et al. 1986, Garland & Else 1987). Less 
is known about individual variation in insect mitochon¬ 
drial activity, but in one study including 19 insect 
species, citrate synthase in flight muscle varied as 
much as three-fold even with only three individuals 
from each species (Alp 1976). Functional correlates of 
variation in mitochondrial enzyme activities have been 
demonstrated in a variety of organisms, most com- 



Volume 57, Number 4 


285 


monly in terms of response to exercise. In a variety of 
mammals, amphibians and reptiles, correlations have 
been seen between mitochondrial activity and maxi¬ 
mum aerobic speed, endurance, and degree of im¬ 
provement after exercise training (Longphre & Gatten 
1994, Garland & Else 1987, Cummings 1979). 

Variables which might be related to mitochondrial 
activity include seasonality and physical characteristics. 
Seasonal variation in mitochondrial enzyme activities 
has been reported in several species, including agamid 
lizards (Garland & Else 1987), iguanid lizards (John- 
Alder 1984), turtles (Olson 1987) and snapper (Majed 
et al. 2002); in all cases the changes are thought to re¬ 
flect seasonal patterns of activity and energetic re¬ 
quirements. In several organisms including lizards and 
fish, mitochondrial enzyme activity is related to body 
size (Garland & Else 1987, Somero & Childress 1980). 
Wing damage may be related to activity levels, as more 
damage is likely to occur during activities such as mat¬ 
ing. Higher wing damage has been seen in male mon- 
archs attempting to mate or drinking dew, as com¬ 
pared to those remaining in clusters (Frey et al. 1998, 
Oberhauser & Frey 1999, Frey et al. 2002). It is possi¬ 
ble that these males represent a more active subset of 
the population, which could be related to metabolic 
differences such as mitochondrial enzyme activity. 

Characterization of mitochondrial activity in overwin¬ 
tering monarchs in central California will provide basic 
information about metabolic variation in this popula¬ 
tion, which may be important for understanding critical 
variables for successful overwintering and reproduction. 
The purposes of this study were to determine the extent 
of variation in succinate dehydrogenase (SDH) activity 
in the male overwintering monarch population and in¬ 
vestigate possible correlations with season or physical 
characteristics. 

Materials and Methods 

This study was carried out at an overwintering site at 
Pismo Beach State Park, Pismo Beach, California. 
Male monarch butterflies were collected monthly at 
approximately 0800 h on the following dates: 20 Octo¬ 
ber, 26 October, 21 November, and 16 December 
2001; 17 January, 19 February and 14 March 2002. In 
all cases this was before the temperature reached the 
flight threshold, so all monarchs were still in clusters 
(Frey et al. 2002). Butterflies were collected from clus¬ 
ters using a standard butterfly net on a long pole, and 
20 males were randomly selected from the net. The 
butterflies were placed in small ziplock bags and 
stored in a styrofoam cooler with icepacks until they 
reached the lab at California Polytechnic State Uni¬ 
versity, approximately fifteen miles away (Frey et al. 


2002). At the lab we measured weight to the nearest 
milligram and forewing length to the nearest millime¬ 
ter from the thorax to the longest extension on the 
forewing. We recorded the number of damaged wings 
containing rips, punctures, or parts missing, and as¬ 
signed each butterfly a wing condition value of 1-4 
based on lack of scales, brightness of color, and dam¬ 
age, with 1 being the worst condition and 4 the best. 
Each butterfly was assigned an identification number 
and placed in a ziplock bag to be stored at -60°C for 
an average of 29 days. There was no relationship be¬ 
tween storage time and enzyme activity (regression 
analysis, R 2 = 4.8%, n = 118)/ 

Succinate dehydrogenase (SDH) activity was mea¬ 
sured in a sample ol thorax muscle tissue from each 
butterfly using the rate of reduction of the artificial elec¬ 
tron acceptor 2,6-dichlorophenolindophenol (DCIP) 
by procedures modified from Singer and Kearney 
(1957). Muscle was obtained by removing all ap¬ 
pendages from the thorax, freezing it in liquid nitro¬ 
gen, and then removing 9-11 mg thorax muscle. The 
tissue sample was weighed and then homogenized 
with a ground-glass homogenizer for 5 minutes on ice 
in 250 pi homogenizing buffer (0.3 M mannitol, 0.02 
M phosphate, pH 7.2). The homogenate was cen¬ 
trifuged at 600 g for 10 minutes at 4°C to remove par¬ 
ticulates and large organelles. Each butterfly tissue 
was assayed in triplicate, using ninety-six well mi¬ 
croplates and a SPECTRAmax Microplate Spec¬ 
trophotometer (Molecular Devices Corp., Sunnyvale 
CA). Reactions were carried out in a total volume of 
200 pi with 30 pi tissue homogenate, 1.5 x 10-4 M 
DCIP, 0.002 M sodium azide and 0.01 M sodium suc¬ 
cinate in assay medium (0.3 M mannitol, 0.02 M phos¬ 
phate, 0.01 M potassium chloride, 0.005 M magne¬ 
sium chloride, pH 7 .2). Reactions were started by 
adding diluted homogenate to the other reagents, and 
a kinetic assay was immediately run at 600 nm every 
15 seconds for 10 minutes. Only the initial linear data 
was used for calculations, typically the first 3 minutes. 
SDH activity in pmoles-min'Tg tissue" 1 was calculated 
from the slope of the initial reaction, the extinction co¬ 
efficient of DCIP (19,100), the measured path length, 
and the volumes and weights used (activity = slope x 
1/e x 1/b x 60sec/min x assay volume/sample homog¬ 
enate vol x total homogenate volume/tissue weight x 
1000 mg/lg). All statistical analyses were performed 
using Minitab (Minitab Inc.). 

Results 

Succinate dehydrogenase (SDH) activity was deter¬ 
mined for a total of 120 butterflies. The overall dis¬ 
tribution of SDH activity for the overwintering season 
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Fig. I. SDH activity in male overwintering monarchs. SDH ac¬ 
tivity is given as flmoles-min^-g tissue -1 . A, Frequency distribution of 
SDH activity for the entire sampled population. B, Boxplot of 
monthly SDH activity: the horizontal line represents the median and 
the ends of the box represent the upper and lower quartiles. 

is shown in Fig. 1A, and ranged from 0.23-3.40 
jLimoles-miir^g tissue -1 , approximately a fifteen-fold 
difference. Ninety percent of the values fall between 
0.5 and 2.7 |Limoles-min -1 -g tissue -1 , a five-fold range. 
Since SDH activity was measured in triplicate samples 
from each butterfly, the standard error associated with 
the mean SDH value for each gives some indication of 
the analytical error; the mean standard error is 0.16, 
approximately eleven percent on average. Given this 
degree of analytical uncertainty, a conservative esti¬ 
mate stills yields at a minimum a four-fold range for 
ninety percent of the enzyme activities and a twelve¬ 
fold range for all values. To illustrate the degree of 
variation in SDH activity, we calculated the coefficient 
of variation (CV = SDx - M00). The overall CV for the 
population over the entire overwintering period was 
forty-seven percent while monthly CVs ranged from 
thirty-six percent in February to fifty-two percent in 
November (Table 1). 

Monthly distributions of SDH activity are shown in 
Fig. IB and Table 1; there was clearly considerable 
variation within each month. Season does have a sig¬ 
nificant association with SDH activity, as seen in a one¬ 


Table 1. Basic statistics for SDH activity, body weight, and wing 
length. SDH activity is given as (J.moles-min _1 -g tissue -1 . Coefficient 
of variation (CV) is the standard deviation divided by the mean, mul¬ 
tiplied by 100. Subscripts (a, b) for SDH activity means are based on 
Tukey s pairwise comparisons; the only significant differences are be¬ 
tween October and November, and November and February. Sam¬ 
ple sizes are n = 20 for each month except for February n = 19. 



Mean 

SD 

Mini¬ 

mum 

Maxi¬ 

mum 

CV(%) 

SDH activity 

October 

1.73 

0.77 

0.57 

3.40 

44 

November 

1.07; 

0.55 

0.23 

2.23 

52 

December 

l-41 ab 

0.58 

0.60 

2.72 

41 

January 

l-<b 

0.73 

0.43 

2.98 

51 

February 

1.80 

0.66 

0.73 

2.82 

36 

March 

l-42 ab 

0.71 

0.49 

2.86 

50 

All butterflies 

1.48 

0.70 

0.23 

3.40 

47 

Body weight (mg) 

October 

599 

80 

454 

692 

13 

November 

536 

68 

490 

705 

13 

December 

529 

73 

395 

745 

14 

January 

539 

68 

434 

670 

13 

February 

557 

65 

430 

660 

12 

March 

478 

87 

330 

650 

18 

All butterflies 

539 

81 

330 

745 

15 

Wing length (mm) 

October 

51.1 

3.0 

45 

56 

6 

November 

50.7 

2.5 

45 

55 

5 

December 

51.9 

2.3 

46 

56 

4 

January 

50.9 

2.1 

47 

55 

4 

February 

48.2 

2.8 

42 

53 

6 

March 

50.1 

2.3 

45 

53 

5 

All butterflies 

50.5 

2.7 

42 

56 

5 


way ANOVA for SDH activity vs. month, F = 3.33, p = 
0.008 (Table 2). Post-hoc pairwise comparisons re¬ 
vealed that the only months with significant differ¬ 
ences are October vs. November and November vs. 
February (Table 1). In terms of seasonal pattern in the 
means, we regressed SDH activity on centered month 
data up to a fifth order polynomial, and only the third 
order polynomial was significant (p = 0.001); however 
this pattern explains very little of the variation (R 2 = 
9.8%). 

Differences in the population at the beginning and 
end of the season may affect the results: in October 
only about 10% of the butterflies have arrived com¬ 
pared to the peak population in December and Janu¬ 
ary, and by March less than 10% of the butterflies are 
left at the overwintering site (D. Frey pers. com.). 
Therefore, we also analyzed the data for the core over¬ 
wintering season when the population is more com¬ 
plete. Consideration of the centered data only from 
November through February reveals a linear pattern 
(p = 0.001), but this still only explains a fraction of the 
individual variation ( R 2 = 13.5%). The estimated in¬ 
crease in SDH activity per month is 0.22 ± 0.06 
jLtmolesmin -1 -g tissue -1 . 






























































Volume 57, Number 4 


287 


Table 2. ANOVA for a general linear model of SDH vs. month 
and all physical variables (body weight, wing length, wing damage, 
and wing condition). Calculations were performed for the full model 
and for month alone, with all physical variables removed. 


Source 

df 

SS 

MS 

F 

P 

Month 

5 

7.4036 

1.4807 

3.33 

0.008 

Physical variables 

4 

1.8438 

0.4610 

1.04 

0.391 

Error 

109 

48.4266 

0.4443 

— 

— 

Total 

118 

57.6740 

— 

- 

— 


We examined two measures of size (body weight 
and wing length) and two measures of wing damage 
(number of damaged wings and overall wing condi¬ 
tion) for possible correlation with SDH activity; size 
data is summarized in Table 1. Mean wet body weight 
was found to decline from October to March, with 
considerable variation apart from the linear decrease 
(simple linear regression analysis, R 2 =10.4%, t = 3.69, 
p = 0.000). Mean wing length fluctuated through the 
overwintering period and declined slightly in February 
and March (simple linear regression analysis, R 2 = 


6.0%, t = 2.74, p = 0.007). The coefficient of variation 
was fifteen percent for body weight and five percent 
for wing length, and the CVs within each month are 
comparable to the overall population. None of the 
physical variables contributed significantly to SDH ac¬ 
tivity (comparison of ANOVA general linear model 
with all variables to model with only month, F = 1.04, 
p = 0.391; Table 2). 

Discussion 

Interindividual variation in succinate dehydro¬ 
genase activity. We found a large degree of interindi¬ 
vidual variation in SDH activity in the male monarch 
population overwintering at Pismo Beach, at least a 
twelve-fold range of enzyme activities. This degree of 
variation is greater than that seen in other physical 
characteristics measured; the coefficient of variation 
for SDH activity is forty-seven percent compared to a 
CV of fifteen percent for body weight and five percent 
for winglength. 

Some of the calculated variation in SDH activity 
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Fig. 2. SDH activity and morphological variables. SDH activity is given as pmoles-min -1 -g tissue -1 . A, SDH activity and butterfly body 
weight. B, SDH activity and butterfly wing length. C, SDH activity and the number of wings damaged. D, SDH activity and wing condition on 
a scale from 1 (worst) to 4 (best). 
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(|imoles*min _1 -g tissue' 1 ) may be an artifact of normal¬ 
izing by tissue weight. Enzyme activities are most 
commonly normalized by wet tissue weight, but dry 
tissue weight, protein content, DNA content, and en¬ 
tire organ weight have also been used (Pelletier et al. 
1995, Cascarano et al. 1978, Spina et al. 1996, Long- 
phre & Gatten 1994). Monarchs are known to vary in 
their state of hydration; in Mexico and southern Cali¬ 
fornia more dehydrated males are seen at the end of 
the overwintering season (Chaplin & Wells 1982, 
Calvert & Lawton 1993), while this pattern has not 
been seen in central California (D. Frey pers. com.). 
Individual differences in hydration status may have in¬ 
troduced some variation into individual SDH activity 
values, but there is no overall correlation between 
lower body weight and higher SDH activity; such a 
correlation might be expected if dehydration led to a 
decrease in body weight but an increase in the number 
of cells per tissue weight. Normalization by protein 
content would address this issue but could introduce 
other problems, since lipid reserves may be depleted 
at the end of the season leading to a breakdown in pro¬ 
tein (Chaplin & Wells 1982, Masters et al. 1988, 
Alonso-Mejia et al. 1997). Future studies will investi¬ 
gate the effects of different normalization approaches 
and possible changes in protein content at the end of 
the overwintering season. 

The degree of interindividual variation in SDH ac¬ 
tivity seen here in male monarchs is within the ranges 
found in other organisms; more importantly, it is large 
enough to have functional consequences since as little 
as a two-fold difference has been correlated with per¬ 
formance differences (Holloszy & Booth 1976). 

Relationship between succinate dehydroge¬ 
nase activity and other variables. Seasonal effects 
are small in comparison to the large variation in indi¬ 
vidual SDH activity levels in butterflies collected on 
the same date. Mean enzyme activity in November is 
significantly lower than in either October or February, 
and there are no significant differences between any 
other months. It is possible that the higher mean SDH 
activity in October could reflect earlier arrival by a 
subset of butterflies with higher mitochondrial activity. 
An increase in mitochondrial activity from November 
to February could be adaptive since the strenuous ex¬ 
ertion of mating occurs late in the season (Tuskes & 
Brower 1978); an increase could also be a conse¬ 
quence of prior flight activity, analogous to the exercise 
training effects seen in other studies. The March mean 
is not significantly different from Februaiy, and the 
downward trend seen could reflect preferential emi¬ 
gration of butterflies with higher mitochondrial ac¬ 
tivity, since only a small fraction of the overwintering 


population remains in March. It is also possible that 
mitochondrial activities in November were lower than 
in other months because of unknown environmental 
factors; preliminary analysis of air temperature data 
did not reveal any correlations with SDH activity or un¬ 
usual occurrences in November (data not shown), but 
careful analysis of potential weather factors has not 
been done. Further investigation will be necessary to 
confirm that the seasonal patterns are seen consistently 
in different years and to determine possible causes. 

None of the physical characteristics examined ap¬ 
pear to be important in determining mitochondrial en¬ 
zyme activity in male overwintering monarchs. Other 
studies have found more linear declines in body 
weight during the overwintering season, but the de¬ 
gree of weight loss depends on the specific overwin¬ 
tering colony examined (Tuskes & Brower 1978, 
Chaplin & Wells 1982, Calvert & Lawton 1993), and it 
is possible that the Pismo Beach site allows better 
maintenance of body weight. The pattern of wing 
lengths seen is similar to that reported by Calvert and 
Lawton (1993) at Mexican overwintering sites except 
that they observed a more dramatic drop in late Feb¬ 
ruary and March after stability for the rest of the sea¬ 
son. They suggest that the decrease at the end of the 
season could be due to larger butterflies leaving the 
colony first. Our observation that mean body weight is 
lower in March but wing length is not, suggests that 
the decreased weight may be due to depleted lipid re¬ 
serves, dehydration, and/or breakdown of protein in 
tissues. Since each of these factors could affect mito¬ 
chondrial enzyme activities, the SDH activity data for 
March are less easily interpreted than for the rest of 
the season. The lack of a correlation with wing damage 
suggests that if there are subsets of the male monarch 
population with different behaviors and corresponding 
degrees of wing damage, mitochondrial activity is not a 
critical determinant nor is it significantly affected by 
behavioral differences which lead to differing wing 
damage. Since we have not examined activities di¬ 
rectly, it is still possible that differing mitochondrial 
enzyme levels do affect activity level, performance or 
mating success. 

Mitochondrial activity could be affected by other 
environmental or physical variables. Age can affect mi¬ 
tochondrial activity and flight performance; changes in 
mitochondrial structure, an increase in mitochondrial 
damage and changes in levels of some metabolic en¬ 
zymes have been reported with aging in the flight mus¬ 
cle of a variety of insect species (Sohal 1976, Ross 
2000). We do not have data on the age of the monarchs 
in this study, but they likely vary by at least a month 
based on emigration data from late summer popula- 
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tions (K. Oberhauser pers. com.). However, we found 
no evidence for an effect of aging on SDH activity 
since the entire population is aging considerably dur¬ 
ing the overwintering period and there was not a uni¬ 
directional seasonal trend. 

Conclusions. The primary conclusion of this study 
is the high degree of interindividual variation in ac¬ 
tivity of a mitochondrial enzyme in the male monarch 
butterflies overwintering at Pismo Beach on the cen¬ 
tral California coast. This variation is not explained by 
seasonality or body size and is not related to the de¬ 
gree of wing damage. Individual mitochondrial en¬ 
zyme activity variation may be partly due to genetic, 
nutritional, and behavioral (especially in regards to 
previous activity levels) differences. Mitochondrial ac¬ 
tivity has potential functional consequences for flight 
performance, and in the case of monarch butterflies, 
potential consequences for survival and reproduction. 
The existence of such substantial individual variation 
suggests that this could be an important factor in indi¬ 
vidual performance and success. Future experiments 
should investigate possible correlations with flight per¬ 
formance and examine mitochondrial activity in fe¬ 
male monarchs, especially with regard to energetically 
demanding reproductive development. In addition, 
further analysis of metabolic parameters at the end of 
the overwintering season would provide valuable infor¬ 
mation about the requirements for successful survival 
and reproduction during the overwintering period. 
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